Evolving Apollonian networks with small-world scale-free topologies 
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We propose two types of evolving networks: evolutionary Apollonian networks (EAN) and gen- 
eral deterministic Apollonian networks (GDAN), established by simple iteration algorithms. We 
investigate the two networks by both simulation and theoretical prediction. Analytical results show 
that both networks follow power-law degree distributions, with distribution exponents continuously 
tuned from 2 to 3. The accurate expression of clustering coefficient is also given for both networks. 
Moreover, the investigation of the average path length of EAN and the diameter of GDAN reveals 
that these two types of networks possess small-world feature. In addition, we study the collective 
synchronization behavior on some limitations of the EAN. 
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I. INTRODUCTION 

In the last few years, complex networks have attracted 
a growing interest from a wide circle of researchers 
[E M IE El : with particular focus on the following three 
properties: power-law degree distribution ;5j, small aver- 
age path length (APL) and high clustering coefficient @. 
The reason for this boom is that complex networks de- 
scribe many systems in nature and society which sharing 
the above-mentioned three characteristics. 

In order to mimic the real-world systems,a wide vari- 
ety of models have been proposed 0, 0, IE 0] • The first 
successful attempt is the Watts and Strogatz's small- 
world network model (WS model) gg, which started 
an avalanche of research on the models of small-world 
networks and the WS model H H EE EH]. An- 
other well-known model is Barabasi and Albert's elegant 
scale- free network model (BA model) which has at- 
tracted an exceptional amount of attention within the 
physics community. In addition to analytic studies of 
the BA model and research of its extensions or modifi- 
cations [l2l Il3l fl4| , many authors have developed a con- 
siderable number of other models and mechanisms that 
may represent processes more realistically taking place 
in real- world networks EE EE E3, EE H3, . 

Until now, modeling complex networks with small-world 
and scale-free properties is still an important issue. 

Recently, based on the well-known Apollonian pack- 
ing, Andrade et al. introduced Apollonian networks |2J| 
which were also proposed by Doye and Massen in |25j 
simultaneously. Apollonian networks belong to a deter- 
ministic growing type of networks, which have drawn 
much attention from scientific communities prl I2H I2H 
H HE HI HI HI H H3. The effects of the Apollo- 
nian networks on several dynamical models have been 
intensively studied, including Ising model and a magnetic 
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model et al HE HE HE] • Doye and Massen adopted 
an extension [2£| of two-dimensional to investigate energy 
landscape networks [39|,|40|. Zhang et al proposed a mini- 
mal iterative algorithm for constructing high dimensional 
networks and studied their structural properties |4lj |. In 
|42| Zhou et al proposed a simple rule that generates 
random two-dimensional Apollonian networks, which are 
generalized by Zhang et al to high dimension [43J. The 
deterministic Apollonian networks (DAN) l24l l25l ItH ] 
and random Apollonian networks (RAN) |42l |43| may 
provide valuable insights into the real-life networks. 

In this paper, first we propose a general scenario for 
constructing evolutionary Apollonian networks (EAN) 
controlled by a parameter q. The EAN can also result 
from Apollonian packing and unifies the DAN and RAN 
to the same framework, i.e., the DAN and RAN are spe- 
cial cases of EAN. Then, we present a general determin- 
istic Apollonian network (GDAN) model governed by a 
single parameter m. Both EAN and GDAN have a power- 
law degree distribution, a very large clustering coefficient 
and a small intervertex separation. The degree exponent 
of EAN and GDAN is changeable between 2 and 3. More- 
over, we introduce an algorithm for the DAN and RAN 
which can realize the construction of our networks. In 
the end, through eigenvalue spectrum of the Laplacian 
matrix, the synchronizability on some limiting cases of 
the EAN is discussed. 



II. BRIEF INTRODUCTION TO 
DETERMINISTIC AND RANDOM 
APOLLONIAN NETWORKS 

We first introduce Apollonian packing (see Fig. 1 
for an example of two dimension), which dates back 
to ancient Greek mathematician Apollonius of Perga. 
The classic two-dimensional Apollonian packing is con- 
structed as follows. Initially three mutually touching 
disks are inscribed inside a circular space which is to be 
ediReip^S. (Jffihifl^rstices °f the initial disks and circle are 
curvilinear triangles to be filled. We denote this initial 



FIG. 1: (Color online) A two-dimensional Apollonian packing 
of disks. 



configuration by generation t = 0. Then in the first gen- 
eration t = 1, four disks are inscribed, each touching all 
the sides of the corresponding curvilinear triangle. The 
process is repeated indefinitely for all the new curvilin- 
ear triangles. In the limit of infinite generations, we ob- 
tain the well-known two-dimensional Apollonian packing. 
The translation from Apollonian packing construction to 
Apollonian network generation is quite straightforward. 
Let the nodes (vertices) of the network correspond to the 
disks and make two nodes connected if the correspond- 
ing disks are tangent 0i US ■ Fig. 2 shows a network 
based on the two-dimensional Apollonian packing. The 
two-dimensional Apollonian network can be easily gener- 
alized to high-dimensions (<i-dimensional, d > 2) [23. El) 
associated with other self-similar packings 0]. The d- 
dimensional Apollonian packings start with d + 1 mu- 
tually touching d-dimensional hyperspheres that is en- 
closed within and touching a larger d-dimensional hyper- 
spheres, with d + 2 curvilinear d-dimensional simplex (d- 
simplex) as their interstices, which are to be filled in suc- 
cessive generations. If each d-hypersphere corresponds to 
a node and nodes are connected if the corresponding d- 
hyperspheres are in contact, then d-dimcnsional Apollo- 
nian networks are established. In every generation of the 
d-dimensional Apollonian packings, if we add only one 
d-hypersphere inside a randomly selected interstice, then 
we get a high dimensional random Apollonian packings, 
based on which high dimensional random Apollonian net- 
works are constructed 0>E3- 



III. THE ITERATIVE ALGORITHMS FOR THE 
NETWORKS 

Before introducing the algorithms we give the follow- 
ing definitions on a graph. The term size refers to the 
number of edges in a graph. The number of nodes in 



FIG. 2: (Color online) Illustration of the two-dimensional 
deterministic Apollonian networks, showing the first two steps 
of iterative process. 



the graph is called its order. When two nodes of a graph 
are connected by an edge, these nodes are said to be ad- 
jacent, and the edge is said to join them. A complete 
graph is a graph in which all nodes are adjacent to one 
another. Thus, in a complete graph, every possible edge 
is present. The complete graph with d graph nodes is 
denoted as (also referred in the literature as d- clique; 
see Two graphs are isomorphic when the nodes of 

one can be relabeled to match the nodes of the other in 
a way that preserves adjacency. Hence all G?-cliques are 
isomorphic to one another. 



A. Evolutionary Apollonian networks 

Now we introduce the evolving Apollonian networks. 
First we give a new packing method for high-dimensional 
(d-dimcnsional, d > 2) Apollonian packings. The initial 
configuration is the same as the deterministic Apollo- 
nian packings. Then in each subsequent generation, each 
d-simplex is filled with probability q. In a special case 
q = 1, it is reduced to the classic deterministic Apollo- 
nian packings, if q approaches but is not equal to 0, it 
coincides with the random Apollonian packings described 
in [4^. l43j . The EAN is derived from this new packing: 
nodes represent d-hyperspheres and edges correspond to 
contact relationship. Fig. 2 shows the network growing 
process for a special case of d = 2 and q = 1. 

In the construction process of the new high- 
dimensional Apollonian packings, for each new d- 
hypersphere, d + 1 new interstices are created. When 
building networks, it is equivalent that for each new 
added node, d + 1 new d-simplices are generated, which 
may create new nodes in the subsequent generations. Ac- 
cording to this, we can introduce a general algorithm to 
create EAN, using which one can write computer pro- 
gram conveniently to simulate the networks and study 
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their properties. 

The d-dimensional EAN after t generations are de- 
noted by A(d,t), d > 2,t > 0. Then at step t, the d- 
dimensional EAN is constructed as follows: For t = 0, 
A(d,0) is a complete graph K4+2 (or (d + 2)-clique). For 
t > 1, A(d, i) is obtained from A(d, £ — 1). For each of 
the existing subgraphs of A(d,t — 1) that is isomorphic 
to a (d + l)-clique and has never generated a node before 
(we call them active (d+ l)-cliques), with probability q, a 
new node is created and connected to all the nodes of this 
subgraph. The growing process is repeated until the net- 
work reaches a desired order. Whe n q = 1. the networks 
are exactly the same as the DAN 0|25|,E|. If q < 1, 
the networks are growing randomly. Especially, as q ap- 
proaches zero and does not equal to zero, the networks 
are reduced to the RAN studied in detail in |42l |43J. 
See for interpretation. 

Next we compute the size and order of EAN. Note 
that the addition of each new node leads to d + 1 new 
(d + l)-cliques and d + 1 new edges. Then, at step 
1, we add expected L v (l) = (d + 2)q new nodes and 
L e (l) = (d + 2)(d + l)q new edges to the graph. Af- 
ter simple calculations, one can obtain that at U(ti > 1) 
the numbers of newly born nodes and edges are L v (ti) = 
q{d+2)(l+dq) t *- 1 and L e (U) = q(d+l)(d+2)(l+dq) t >- 1 , 
respectively. Thus the average number of total nodes N t 
and edges E t present at step t is 

t 

N t = {d + 2) + J2 L v(ti) 

ti=i 

= ( d W l+dq) \ +d - 1 a) 



and 



respectively. So for large t, The average degree fc t = 
is approximately 2(d + 1). Moreover, when d = 2, we 
have i?t = 37V t — 6. Thus in this case, all networks are 
maximal planar networks (or graphs) |42 |. 



the existing subgraphs of G(d, t — 1) that is isomorphic to 
a (d+ l)-clique and created at step t — 1, m new vertices 
are created, and each is connected to all the vertices of 
this subgraph. In the limit of infinite generations we ob- 
tain GDAN. When m = 1, the networks are reduced to 

dan mmim. 

Let L v (t), L e (t) and Ku+x),t De the numbers of ver- 
tices, edges and (d+ l)-cliques created at step t, respec- 
tively. Because the addition of each new vertex leads 
to d + 1 new (d + l)-cliques and d + 1 new edges, we 
have L v (t) = mify+i),*-]., L e (t) = (d + 1)1^ (t), and 
-^(d+i).t = (d + l)itj(^). Thus one can easily obtain 
if (d+1) ; t = m(d + 1)^+!).^! = (d + 2)[m(d + 1)]* 
(t > 0), L„(t) = (d + 2)m*(d + l)*" 1 (t > 0) and 
L e (t) = (d + 2)m t (d+l)* (t > 0). From these results, we 
can easily compute the size and order of the networks. 
The total number of vertices Nt and edges E t present at 
step t is 



N t = ^n v {U) 

ti=0 

m(d + 2)[m t (d+l) t 
m(d+l)-l 



d + 2 



(3) 



and 



(4) 



E t — n e (ti) 

ti=0 

_ m(d + 2)(d + l)[m*(d + 1)* - 1] (d + 2)(d + 1) 
~ m(d+ 1) - 1 + 2 : 

respectively. So for large large, The average degree k t 
approaches 2(d + 1). 



IV. RELEVANT CHARACTERISTICS OF THE 
NETWORKS 

In the following we will study the topology properties 
of EAN and GDAN, in terms of the degree distribution, 
clustering coefficient, average path length and diameter. 



A. Degree distribution 



B. General deterministic Apollonian networks 

According to the construction process of d-dimensional 
Apollonian packings, in Ref. jjjj, a generation algorithm 
for d-dimensional Apollonian networks was proposed. 
Here we generalize the algorithm to establish general de- 
terministic Apollonian networks (GDAN). The network, 
denoted by G(d, t) after t generations with d > 2 and 
t > 0, is constructed in an iterative way. For t = 0, 
G(d, 0) is a complete graph Kd+2 (or (d + 2)-clique). For 
t > 1, G(d,t) is obtained from G(d,t — 1). For each of 



The degree distribution is one of the most important 
statistical characteristics of a network. We will discuss 
the degree distribution of EAN and GDAN in detail. 

1. Evolutionary Apollonian networks 

When a new node i is added to the graph at step U , it 
has degree d + 1 and forms d + 1 active (d + l)-cliques. 
Let L c (i,t) be the number of active (d+l)-cliques at step 
t that will possibly created new nodes connected to the 
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node i at step t+1. Then at step U, L c (i,ti) = d + 1. 
At step ti + 1, there arc (d + l)q new nodes which forms 
(d+l)qd new active (d+l)-cliques containing i, and there 
are (d+l)q active (d+l)-cliques of i are deactivated at the 
same time. Thus L c (i, ti + 1) = (d + 1)[1 + (d — l)q]. It is 
not difficult to find the following relation: L c (i,t + 1) = 
L c (i,t)[l + (d — l)q]. Since L c (i,U) = d + 1, we have 
L c (i,t) = (d + l)[l + (d-l)g]*-*«. Then the degree h(t) 
of node z at time t is 



jfei(t) = d + l + q^2L c (i,T) 

■ ^( iVg^ , (5) 

Since the degree of each node has been obtained ex- 
plicitly as in Eq. J5J, we can get the degree distri- 
bution via its cumulative distribution, i.e., P C um{k) = 
J2k'>k N ( k '^)/ N t ~ fcl " 7 : where AT(fc',t) denotes the 
number of nodes with degree k' . The detailed analysis is 
given as follows. For a degree k 




fc= (d+1) 



[1 + (d - l)q] t - m + d - 2 



1 



FIG. 3: (Color online) The cumulative degree distribution 
Pcum(k) at various q values for dimension d = 2. The circles 
(a), squares (b), stars (c), and triangles (d) denote the sim- 
ulation results for networks with different evolutionary steps 
t = 1110, t = 23, t = 14, and t = 12, respectively. The 
four straight lines are the theoretical results of -y(d, q) as pro- 
vided by equation (6). All data are from the average of ten 
independent simulations. 



there are L v (m) = q(d + 2)(1 + dg) m_1 nodes with this 
exact degree, all of which were born at step m. All nodes 
born at time m or earlier have this or a higher degree. 
So we have 



E 

k'>k 



N(k',t) = (d + 2)+Vi t (s) 



i 



(d + 2 



(l + dq)" 



d-l 



i.e., q = 1 and q — > (but q ^ 0), the evolutionary 
Apollonian network reduces to the deterministic Apollo- 
nian networks [24l I25L |4~H and their stochastic variants 
^jEHi respectively. Fig. (0 shows, on a logarithmic 
scale, the scaling behavior of the cumulative degree dis- 
tribution P CU m{k) for different values of q in the case of 
d = 2. Simulation results agree very well with the ana- 
lytical ones. 



As the total number of nodes at step t is given in Eq. (JTJ 
we have 



(d+1) 



[1 + (d- l).?]'-" 1 +d- 2 



d-l 
(d-l- 



1-7 



(d + 2) 

Therefore, for large t we obtain 

[[l + (d -I)?]*""*] 1-7 = (1 



(l+dq)'+d-l 



•dg) r 



and 



1 



ln(l + dq) 
ln[l + (d- l)q]' 



(6) 



Thus, the degree exponent 7 is a continuous function 
of d and g, and belongs to the interval (2,3]. For any 
fixed d, as q decrease from 1 to 0, 7 increases from 



1 + MiM to 2 



In the case d 



1 ■ . ..... >..■"< -2,7 can be 

tunable between 2.58496 and 3. In the two limitations, 



2. General deterministic Apollonian networks 

Considering a vertex i added to the networks at step ti. 
Let if (d+i) (i, t) be the number of newly-created (d + 1)- 
cliques at step t containing vertex i. These new cliques 
will create new vertices connected to the vertex i at step 
t+1. At step U, K(d+i)(i,U) = d+ 1. From the iterative 
process, one can see that each new neighbor of i generates 
d new (d + l)-cliques with i belonging to them. Then it 
is not difficult to find following relations: 



Afcj(t) = h{t) - h{t - 1) = mK (d+1) (i,t - 1) 



and 



K (d+1) {i,t) = dmK {d+1) {i,t-l) = {d+l){dm)^ u 

Then the degree of vertex i becomes 

t-i 

h{t) = kj(tj) + m K( d+1 )(i,T) 

r=t< 

m(rf+l)[H) M '-l] + (l 2 -l 



(7) 
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For a degree k 

m(d+ l)[(md) t -' p - 1] + d 2 - 1 
k= —1 ' 

there are n v (p) = (d + 2)m p (d + l^" 1 vertices with this 
exact degree. Also, we have 

k'>k s=0 

m(d + 2)[mP(d+ l) p - 1} 



m(d+l) - 1 



d+ 2. 



(8) 



From the definition of cumulative degree distribution, we 
have 



i(d + l)[{mdy~P - 1] + d 2 - 1 



d-1 



1-7 



m(d+2)[m p (d+l) p -l] , , , 2 
m(d ••!) • 1 
~~ m(rf+2)[m'(d+l) t -l] I ^ I o ' 
m(d+l)— 1 

When i is large enough, one can obtain 

[{mdy-P} 1 " 1 = [m(d+l)f- t 

and 



7 w 1 + 



In[m(ef+1)] 
ln(mci) 



(9) 



For m = 1, Eq. © recovers the results previously ob- 
tained in Refs. EllHEl- 



B. Clustering coefficient 

The clustering coefficient jfj Cj of node i is defined as 
the ratio between the number of edges e, that actually 
exist among the fcj neighbors of node i and its maximum 
possible value, — 1)/2, i.e., C; = 2ej/[fcj(fcj — 1)]. The 
clustering coefficient of the whole network is the average 
of C i s over all nodes in the network. 

For both EAN and GDAN, the analytical expression of 
clustering coefficient C(k) for a single node with degree 
k can be derived exactly. When a node is created it is 
connected to all the nodes of a (d + l)-clique, in which 
nodes are completely interconnected. So its degree and 
clustering coefficient are d+1 and 1, respectively. In the 
following steps, if its degree increases one by a newly cre- 
ated node connecting to it, then there must be d existing 
neighbors of it attaching to the new node at the same 
time. Thus for a node of degree fc, we have 



d{d+l) 



C(k) 



d(k-d-l) 



fc(fc-i) 

2 



2d(k- 
fc(fc-l) ' 



(10) 



which depends on degree k and dimension d. For k 
d, the C(k) is inversely proportional to degree k. The 
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FIG. 4: (Color online) The clustering coefficient of the whole 
network as a function of q and d. Results are averaged over 
ten network realizations for each datum. 



scaling C(k) ~ k 1 has been found for some network 
models [SlJllJIlJsgj^JMIillilEgl, including DAN 
and RAN |24ll25l l4ll 142 . 14^, an d has also been observed 
in several real-life networks [3(| ■ 

Using Eq. 1)10(1 . we can obtain the clustering Ct of the 
networks at step t: 



1 « 2d(D r - ±t±)L v (r) 



(ii) 



where the sum runs over all the nodes and D r is the 
degree of the nodes created at step r, which is given by 
Eq. @or0. 



1. Evolutionary Apollonian networks 

For EAN, it can be easily proved that for any fixed 
d, Ct increases with g, and that for arbitrary fixed q, 
Ct increases with d. Exactly analytical computation 
shows: in the case d = 2, when q increases from to 
1, C grows from 0.7366 to 0.8284, with a special value 
C t = 0.7934 for q = 0.5. Likewise, in the case d — 3, C 
increases from 0.8021 to 0.8852 when q increases from 
to 1, especially C t — 0.8585 for q = 0.5 (see also 
[H> OH [3 EH)- Therefore, the evolutionary networks 
are highly clustered. Fig. shows the clustering coeffi- 
cient of the network as a function of q for d = 2 and d = 3, 
respectively, which is in accordance with our above con- 
clusions. From © and Figs. 10}, one can see that both 
degree exponent 7 and clustering coefficient Ct depend 
on the parameter q. The mechanism resulting in this 
relation deserves further study. The fact that a biased 
choice of the cliques at each iteration may be a possible 
explanation, see Ref. |47| . 
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FIG. 5: (Color online) The dependence relation of C on d 
and m. 



2. General deterministic Apollonian networks 

For GDAN, in the infinite network order limit (N t — + 
oo), Eq. (|llfl converges to a nonzero value. When d = 2, 
for m = 1, 2 and 3, C is equal to 0.8284, 0.8602 and 
0.8972, respectively. When to = 2, for d = 2, 3 and 4, 
C are 0.8602, 0.9017 and 0.9244, respectively. Therefore, 
the clustering coefficient of GDAN is very high. More- 
over, similarly to the degree exponent 7, clustering coef- 
ficient C is determined by both d and to. Fig. \5\ shows 
the dependence of C on d and to. From Fig. [3] (a) and 
(b), one can see that for any fixed to, C increases with 
d. But the dependence relation of C on to (see Fig. [5] 
(b)) is more complex: (i) when m < 2 and d < 6, for the 
same d, C increases with to; (ii) when m < 2 and rf > 6, 
for the same d, C decreases with m; (iii) when to > 3, 
for arbitrary fixed d, C increases with to. Further effort 
should be paid for this complicated relation. 



approximately as: 

a{N + 1) = a(N) + A + (A - d - l)£(N - d). (13) 
It is clear 

(N — d- 1)£(N - d) = — tt— — r-^ < —77-^- (14) 



N-d 



N 



From Eqs. (|13fl and l|14|) . we can provide an upper bound 
for the variation of cr( A) as 



da(N) =N+ 2o-(N) 



dN 



N ' 



which leads to 



a(N) = N 2 {\nN + P), 



(15) 



(16) 



where (3 is a constant. As o~{N) ~ A 2 In TV, we have 
£(A) - In TV. 

It should be emphasized that as Eq. (|15fl has been de- 
duced from an inequality, then £(N) increases at most 
as In A with A. Here we only give an upper bound for 
APL, which increases slightly slower than In A. Thus, 
our model exhibits the small-world property. Especially, 
in the case of q = 1, we can compute exactly the diameter 
of the network, which is the maximum distance between 
all node pairs of a graph. Previously reported analytical 
result has shown that the diameter grows logarithmically 
with the order of the network |4l|. In Fig. ©, we report 
the simulation results on the APL of networks for differ- 
ent q and d. From Fig. JfJJ), one can see that for fixed d, 
APL decreases with increasing q; and for fixed q, APL 
is a decreasing function of d. When network order N is 
small, APL is a linear function of InN; while A becomes 
large, APL increases slightly slower than InN. 



D. Diameter for general deterministic Apollonian 
networks 



C. Average path length for evolutionary 
Apollonian networks 

We label the nodes by their creation times, v = 
1, 2, 3, . . . , A - 1, N. Using i{N) to represent the APL 
of the our networks with order A, then we have: £(N) = 

n(n N \) i where a(N) = Y,i<i<j<N d hi is the total dis " 
tance, in which di j is the smallest distance between node 
i with node j. Now we study the APL of the present 
model by using the approach similar to that in 0> 03 • 
Since the distances between existing node pairs will not 
be affected by the addition of new nodes, thus we have: 

N 

a{N + l)=a{N)+J2di,N+i. (12) 
Like in the analysis of 0> 03 j Eq. (|12|) can be rewritten 



The diameter of a network characterizes the maximum 
communication delay in the network and is defined as 
the longest shortest path between all pairs of vertices. 
In what follows, the notations \x~\ and [^J represent 
the integers obtained by rounding x to the nearest in- 
tegers towards infinity and minus infinity, respectively. 
Now we compute the diameter of A(d,t), denoted by 
diam(A(d,t)) for d > 2: 

Step 0. The diameter is 1. 

Steps 1 to |~|]. In this case, the diameter is 2, since 
any new vertex is by construction connected to a (d+1)- 
clique, and since any (d + l)-clique during those steps 
contains at least | + 2 (d even) or + 1 (d odd) vertices 
from the initial (d + 2)-clique A(d, 0) obtained after step 
0. Hence, any two newly added vertices u and v will 
be connected respectively to sets S u and S v , with S u C 
V(A(d,0)) and S v C V(A(d,0)), where V(A(d,0)) is the 
vertex set of A(d, 0); however, since \S U \ > | + 2 (d even) 
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FIG. 6: Semilogarithmic graph of the APL vs the network 
order N. Besides N, APL depends on d and q. Each data is 
obtained by ten independent network realizations. The lines 
are linear functions of In N. 

and \S V \ > ^- + 1 (d odd), where IS"! denotes the number 
of elements in set S, we conclude that S u H S v 7^ 0, and 
thus the diameter is 2. 

Steps [|] +1 to d+1. In any of those steps, some newly 
added vertices might not share a neighbor in the original 
(d + 2)-clique A(d, 0) obtained after step 0; however, any 
newly added vertex is connected to at least one vertex of 
the initial clique A(d, 0). Thus, the diameter equals to 3. 

Further steps. Clearly, at any step t > d + 2, the di- 
ameter always lies between a pair of vertices that have 
just been created at this step. We will call the newly 
created vertices "outer" vertices. At any step t > d + 2, 
we note that an outer vertex cannot be connected with 
two or more vertices that were created during the same 
step < t' < t — 1. Moreover, by construction no two 
vertices that were created during a given step are neigh- 
bors, thus they cannot be part of the same (d+ l)-clique. 
Thus, for any step t > d + 2, some outer vertices are con- 
nected with vertices that appeared at pairwise different 
steps. Thus, there exists an outer vertex v t created at 
step t, which is connected to vertices w-s, 1 < i < t — 1, 
all of which are pairwise distinct. We conclude that v t 
is necessarily connected to a vertex that was created at 
a step to < t — d — I. If we repeat this argument, then 
we obtain an upper bound on the distance from v t to 
the initial clique A(d, 0). Let t — a(d + 1) + (3, where 
1 < (3 < d+1. Then, we see that v t is at distance at most 
a + 1 from a vertex in A(d, 0). Hence any two vertices 
v t and w t in A(d, t) lie at distance at most 2{a + 1) + 1; 
however, depending on /3, this distance can be reduced 
by 1, since when fj < \^~\ , we know that two vertices cre- 
ated at step (5 share at least a neighbor in A(d, 0). Thus, 
when 1 < p < [|] , diam(A(d, t)) < 2(a + 1), while when 
I"! 1 + 1 < p < d + 1, diam{A(d, 1)) < 2(a + 1) + 1. One 
can see that these distance bounds can be reached by 



pairs of outer vertices created at step t. More precisely, 
those two vertices Vt and Wt share the property that they 
are connected to d vertices that appeared respectively at 
steps t — 1, t — 2, . . . t — d — 1. 

Based on the above arguments, one can easily see that 
for t > d+2, the diameter increases by 2 every d+1 
steps. More precisely, we have the following result, for 
any d > 1 and t > 1 (when t = 0, the diameter is clearly 
equal to 1): 

diam(A(d,t)) = 2([tzl\ +l) + /(d,t) 

where f(d, t) = if t - L^jJ {d+1) < [f ] , and 1 other- 
wise. 

In the limit of large t, diam(A(d,t)) ~ ^j-, while 
N t ~ [m{d + 1)]*, thus the diameter is small and scales 
logarithmically with the network order. 

V. SYNCHRONIZATION ON SOME LIMITING 
CASES FOR EVOLUTIONARY APOLLONIAN 
NETWORKS 

The ultimate goal of the study of network struc- 
ture is to study and understand the workings of sys- 
tems built upon those networks^ 0, El El Recently, 
some researchers have focused on the analysis of func- 
tional or dynamical aspects of processes occurring on 
networks. One particular issue attracting much atten- 
tion is the synchronizability of oscillator coupling net- 
works 48] . Synchronization is observed in diverse natural 
and man-made systems and is directly related to many 
specific problems in a variety of different disciplines. It 
has found practical applications in many fields including 
communications, optics, neural networks and geophysics 
IU HI EH!! Hail. After studying the relevant char- 
acteristics of network structure, which is described in the 
last sections, we will study the synchronization behavior 
on the networks. 

We follow the general framework proposed in [5^, |5^| , 
where a criterion based on spectral techniques was estab- 
lished to determine the stability of synchronized states on 
networks. Consider a network of N identical dynamical 
systems with linearly and symmetric coupling between 
oscillators. The set of equations of motion for the sys- 
tem are 

N 

X S ; =F(x i )+<7^GyH(x i ), (17) 

3=1 

where = F(x^) governs the dynamics of each individ- 
ual node, H(xj) is the output function and a the cou- 
pling strength, and Gij is the Laplacian matrix, defined 
by Ga = ki if the degree of node i is ki, = — 1 if 
nodes i and j are connected, and Gij = otherwise. 

Since matrix G is positive semidefinite and each rows 
of it has zero sum, all eigenvalues of G are real and non- 
negative and the smallest one is always equal to zero. 
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FIG. 7: (Color online) The eigenratio R as a function of 
network order N. All quantities for random networks are 
averaged over 50 realizations. 



We order the eigenvalues as = Ai < A2 < • • • < A at. 
Then one can use the ratio of the maximum eigenvalue 
A at to the smallest nonzero one A? to measure the syn- 
chronizability of the network pSL l56| . If the eigenratio 
R = Xn/\2 satisfies R < 012/ a±, we say the network is 
synchronizable. Here the eigenratio R depends on the 
the network topology, while a.%ja,\ depends exclusively 
on the dynamics of individual oscillator and the output 
function. Ratio R = represents the synchroniz- 

ability of the network: the larger the ratio, the more dif- 
ficult it is to synchronize the oscillators, and vice versa. 

After reducing the issue of synchronizability to finding 
eigenvalues of the Laplacian matrix G, we now investigate 
the synchronization of our networks. Here we only study 
two limiting cases: q = 1 |2jJ |4l| and q — >• (but 
1 7^ 0) |42l |43| . The eigenratio R of different networks is 
obtained numerically for different d and q, as exhibited 
in Fig. J7J). One can see that for the same d and N, 
R of the deterministic networks is smaller than that of 
their corresponding random versions, which implies that 
the synchronizability of the former is better. Even in the 
case of different d, the deterministic networks (d = 2) 
are easier to synchronize than those random networks 
(d = 3) with the same order but different average node 
degree. 

Why coupling systems on the two class of networks 
exhibit very different synchronizability? Previously re- 
ported results have indicated that underlying network 
structures play significant roles in the synchronizability of 
coupled oscillators. However, the key structural feature 
that determines the collective synchronization behavior 
remains unclear. Many works have discussed this issue. 
Some authors believe that shorter APL tends to enhance 
synchronization [5ol IBH l58j. In contrast, Nishikawa et. 
al. reported that synchronizability is suppressed as the 
degree distribution becomes more heterogeneous, even 
for shorter APL [5|J. Also, some investigations showed 



that R decreases whenever the betweenness heterogeneity 
decreases l60l. while an opposite conclusion was claimed 
in |6l|]. In [62j. the authors asserted that larger average 
node degree corresponds to better synchroizability. All 
these may rationally explain the relations between syn- 
chroizability and network structure in some cases, but 
do not well account for the synchroizability on the evolu- 
tionary Apollonian networks. More recently, some other 
authors have presented that structure and distribution of 
hubs is the key in what refers to enhance synchroizabil- 
ity 63, 64], which may be a possible explanation for the 
better synchroizability of deterministic networks when 
compared to random ones. But we speculate that the 
synchroizability on the evolutionary Apollonian networks 
is not determined by a single structure property, but by 
the combination of APL, heterogeneity of degree distri- 
bution, betweenness centrality, modularity, mean node 
degree, and so on |63| , which need further deep research. 



VI. CONCLUSION AND DISCUSSION 

In summary, on the basis of Apollonian packings, we 
have proposed and studied two kinds of evolving net- 
works: evolutionary Apollonian networks (EAN) and 
general deterministic Apollonian networks (GDAN). Ac- 
cording to the network construction processes we have 
presented two algorithms to generate the networks, based 
on which we have obtained the analytical and numeri- 
cal results for degree distribution and clustering coeffi- 
cient, as well as the average path length, which agree 
well with large amount of real observations. The degree 
exponents can be adjusted continuously between 2 and 3, 
and the clustering coefficient is very large. Moreover, we 
have studied the synchronization of some limiting cases 
of the EAN and found that the stochastic networks are 
more hardly synchronized than their deterministic coun- 
terparts. 

Because of their three important properties: power- 
law degree distribution, small intervertex separation and 
large clustering coefficient, the proposed networks pos- 
sess good structural features in accordance with a variety 
of real-life networks. For the special case of d = 2, the 
networks are maximal planar graphs. This may be help- 
ful for designing printed circuits. Moreover, our networks 
consist of complete graphs, which has been observed in 
variety of the real- world networks, such as movie actor 
collaboration networks, scientific collaboration networks 
and networks of company directors 0, 0> 0> • 
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